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The ionization structure of early-B supergiant winds 
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Abstract. We present empirically determined ionization conditions for the winds of 106 luminous BO to B5 stars observed by 
lUE. The UV wind lines are modelled to extract products of mass-loss rates times ionization fractions (M^,(w), where w = 
v/v„) for Nv, Civ, Si IV, Sim, Aim and Cii. We describe the general behaviour of the M qi(w) and their ratios, demonstrating 
that the wind ionization increases with distance from the star, contrary to recent findings for O star winds. Using empirical 
mass-loss rates (from Ha observations) and model prescriptions, we derive mean qj(w) values integrated over the wind, (^,>. 
These (g,) are quite small, never exceeding 15% for Al iii or 2% for Si iv. This is surprising, since the (qi) for these ions clearly 
peak within the observed spectral range. We conclude that the low (^,> arise because the (M qj) are underestimated by the wind 
models, which assume that the outflows are smooth when they are, in fact, highly structured. 
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1. Introduction 

Luminous massive stars are key constituents of galaxies since 
they provide prolific ionizing radiation and an input of energy 
and matter to the surrounding environment through strong ra- 
diatively driven winds. More specifically, early-type B super- 
giants (Teff ~ 12 to 30kK) are the most numerous hot lumi- 
nous stars. They make substantial contributions to the inte- 
grated light of starbursts and provide diagnostics of star for- 
mation at low and high redshifts (e.g. de Mello et al. 2000). B 
supergiants are also key markers in the determinations of extra- 
galactic distances based on the Wind-Momentum-Luminosity 
(WML) relation (e.g. Kudritzki & Przybilla 2003). 

With this in mind, it is clearly important to understand the 
mass-loss process in B supergiants and derive reliable parame- 
ters for mass-loss rates, M, ionization mixture and wind struc- 
ture. In this Letter we focus on the empirical ionization struc- 
ture, as part of a wider comprehensive survey of B supergiants 
which will be published in a forthcoming paper (Searle et al., 
in preparation). Line-synthesis modelling has been employed 
to derive the wind properties of a large sample of B0-B5 stars, 
for ions with resonance lines available in high-resolution lUE 
data. Our primary aim here is to examine the properties of the 
ionization state of the winds, determine whether specific ions 
are near the dominant stage of ionization, and examine how 
the properties change with spectral type and distance from the 
star. Our empirical results provide critical constraints for more 
detailed wind modelling. We highlight the unexpectedly low 
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values that are apparent for the ionization fractions and the po- 
tential role of wind clumping. 

2. SEI modelling of lUE spectroscopic data 

Our study of Galactic B supergiant winds is based on high- 
resolution (A/AA ~ 10'^) WE SWP (1150 < A < 1900A) 
spectra of 106 stars. The individual spectra were binned on a 
wavelength grid with regular sampling of 0. 1 A, and the typical 
signal-to-noise in the continuum is ~ 20. The UV spectra of 
B supergiants contain wind lines from a wide range of ioniza- 
tion states, including Cii ^.11334.53, 1335.71, Civ .1/11548.20, 
1550.77, Nv/I^1238.82, 1242.80, Al iiU^ 1854.72, 1862.79, 
Si III /l/l 1206.50 and Si iv AAU92>.16, 1402.77. 

To extract physical parameters from the wind profiles, we 
used the methods described in detail by Massa et al. (2003). 
They employ a modified version of the 'Sobolev with exact in- 
tegration' (SEI) code (Lamers et al. 1987). Once the velocity 
field of the wind is determined, this approach provides reliable 
radial optical depths for unsaturated wind lines as a function 
of normalized velocity, t,o^(w), where w - v/vco and Voo is the 
terminal velocity of the wind. The fits to individual wind lines 
require parameters for the terminal velocity, 'yS-type' velocity 
law, turbulent velocity, and t„^(w) (specified as a set of 21 in- 
dependent velocity bins that are adjusted by a non-linear least 
squares procedure). Photospheric spectra were supplied using 
the lUE line profiles of the non-supergiant B stars listed by 
Prinja et al. (2002). The T,ad{w) were converted into a product 
of M times ionization fraction, M qi{w) (see, e.g., Massa et al. 
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Fig. 1. Left panels: (Mqi) values (integrated over 0.2 < v/voo < 0.8) versus T^b for the ions studied. Right panel: (qt) versus Jeff 
for the same ions. Solid points are measurements from unsaturated lines of stars with logCL/L©) > 5, crosses are for similar lines 
from stars with logCL/L©) < 5. Upward arrowheads indicate lower limits determined from saturated lines ((T„,f/) > 5 for Si^^ and 
C^^, 6 for C^^ and 7 for N'*^, Si^^ and Al^^), and downward arrowheads indicate upper limits from non-detections (Tradbv) < 0.1 
for all w). 



2003). To isolate qi(w) for each ion requires an estimate of M. 
We discuss our approach to this problem in § 3.1. 

3. Ionization conditions 

Figure 1 (left panels) shows mean {M q,) values over 0.2 < 
v/Vcx, < 0.8, for the ions analysed versus T^.^, and (right panels) 
the mean {qi) values for the same ions (see § 3.1). We adopted 
the Humphreys & McElroy (1984) spectral type - Teff calibra- 
tion. Recent detailed NLTE model atmosphere analyses (e.g. 
Crowther et al. 2002, Bianchi & Garcia 2002, Herrero et al. 
2002) have revised the T^s scale for O- stars downward by ~ 



10 - 20%. There are also indications that B star temperatures 
should be reduced as well (e.g. Evans et al. 2004, Searle et 
al. 2005). However, while the revised Teff of a BO la may be 
~ 2000 K lower than those used in Fig. 1, the differences are 
expected to be 5 5% at Bl and later (Searle et al. 2005). 

We first consider the {Mqd- These behave as expected: at 
fixed Teff, all of the (Mqi) tend to be larger for more lumi- 
nous stars, with some overlap, and the presence of an ion at a 
given Teff is related to the energies required to produce and de- 
stroy it. Specifically, {Mq(N^'^)} and {Mq(C^'^)) increase over 
the entire Tes range. (M^(Si^'^)) appears to peak (or at least 
plateau) at ~ 20kK (with considerable spread at each T^s). The 
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Fig. 2. Ionization fraction ratios of Si /N (upper panel) and 
^2+^(^3+ (loyygj. panel) plotted as a function of normalised ve- 
locity (w = v/Voo). Only the cases of unsaturated line profiles 
are shown. 



peak is more distinct in the less luminous stars, where satu- 
ration is not a problem. {M q(Si^^)} (which can be measured 
in fewer stars) peaks at ~ 20kK. {Mq(Al^'^)} peaks at 18kK 
with a value ^ 3 x 10"^ Mq yr"', the largest of all ions. Aim 
is a useful ion for mass-loss studies of cooler B stars, since 
it rarely saturates (due to its relatively low abundance) and is 
unafi^ected by nuclear processing. Note that {M q(Al^'^)} peaks 
at a lower T^g than {M q{Si^^)}, even though a higher energy 
is required to produce it. This occurs because it is destroyed 
at a lower energy than Si^"^. {Mq(C^)) peaks at 18kK and then 
abruptly drops to zero - independently of the stellar luminosity. 
This unique behaviour may be due to the near coincidence of 
the ionization potentials of C+ (24.38 eV) and He (24.48 eV). 
When He ionizes in the stellar photosphere, the stellar radiation 
field in the He"^ continuum (228 <A< 506 A) increases dramat- 
ically. Indeed, Kurucz (1991) logg = 3 models with T^s =17 
and 19 kK show that the integrated flux in this region increases 
by a factor of 30. This enormous increase could convert all of 
the C"^ in the wind to C^"^. However, these conjectures must be 
verified by detailed modelUng. 

We now consider M q ratios, which are independent of M. 
We begin with ratios of the (M qi), which provide insight into 
the how the wind ionization responds to T^s- The overall trend 
is that the empirical ionization ratios of N^"^/Si^"^, C^'^/Si^"^ 
and Si^'^/Al^'^ increase as a function of T^s, indicating that the 
winds of hotter B supergiants are more highly ionized. Some of 
the ratios straddle the 'bistability jump' at ~21kK, where there 
is a ramped increase in the ratio of terminal velocity to escape 
velocity (e.g. Lamers et al.l995, Vink et al. 2000). However, 
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Fig. 3. Mean ionization fraction of Si (filled) and Al (open), 
based on the Ha mass-loss rates of Kudrtitzi et al. (1999; K99). 

the ratios of C^^/Si^^, and Si^^/Si^"^ do not suggest a dramatic 
change in the ionization of the winds across the 'jump', though 
the former has a steeper drop below ~ 20kK. 

FinaUy, we examine how the ionization of B supergiant 
winds varies with distance from the star by plotting M qi(w) 
ratios for Si-'^N'*^ and Al^^/C^^ against w (= v/voo) in Fig. 2 
(for unsaturated cases only). Broadly, the ratios decrease over 
0.2 < w < 0.9 (typically 1.2 < r//?* < 10), implying an in- 
crease in the ionization state with velocity. This is the oppo- 
site of the trend found by Massa et al. (2003) in LMC O stars, 
and is more in agreement with the expectations of the optically 
thin nebular approximation (e.g., Cassinelli & Olson 1979). It 
could indicate a fundamental difference between the winds of 
O and B-type stars. We intend to compare these trends to de- 
tailed model predictions in a future pubUcation. 

3.1. Mean ion fractions and M 

Conversion of our (M qi) measurements into {qi), requires M 
values. We consider two options. First we adopt empirical 
mass-loss rates determined from a NLTE model atmosphere 
analysis of Hor profiles by Kudritzki et al. (1999). This accounts 
for only 14 stars in our sample, but it allows us to compare 
our results to previous work. The resulting ion fractions for 
Aim and Siiv are shown in Fig. 3, where we have adjusted 
the Kudritzki et al. Ms so that the adopted radii and tempera- 
tures are consistent with those used in Fig. 1 (the adjustments 
are typically less than a factor of 2, and not systematic). We see 
that (^(AP+)) peaks at ~ 10-15%, which is comparable to the 
maximum value of {q{^*^)) found by Massa et al. (2003, 2004) 
for O stars mass-loss rates based on the far-UV P v wind lines 
and either empirical or theoretical Ms. The (^(Si^^)) peak is 
poorly sampled by the data and dominated by low luminosity 
stars, since the more luminous stars have saturated Si iv wind 
Unes. 

Second, since the vast majority of our program stars do not 
have empirically determined mass-loss rates, we use the Vink 
et al. (2000, 2001) model prescriptions between stellar param- 
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eters and M. Note that changes of Teff or L/Lq resulting from a 
± one spectral or luminosity bin error, change these predictions 
by less than a factor of 2, and that these values are far larger 
than differences between our adopted parameters and the ones 
determined by Kudritzki et al. (1999) for stars in conmion. The 
right side of Fig. 1 shows the resulting {q;} as a function of Teff. 
Overall, these follow the same trends as the {M qi). While it ap- 
pears that the less luminous stars have larger (g,), this is most 
likely a result of line saturation in the more luminous stars. 

We now consider the four ions which peak in the B star 
range, to determine if any of them are dominant. We find that 
(<7(Si"'^)) never exceeds 1.9%, and that the mean of all unsat- 
urated cases near its peak is 0.48%. These quantities for the 
other ions are: 1.2 and 0.14% (Si^^), 6.9 and 1.7% (Al^^) and 
0.34 and 0.08% (C+). The {q{A\^*)) values are based primar- 
ily on high luminosity stars, and these exceeds 2.5% in only 
one case - somewhat smaller than Massa et al. found for P v 
in the O stars. The (<7(Si^^)) values are particularly surprising, 
since models predict it should be ~ 1 for 18 < Jeff < 20kK. 
However, due to the saturation of the Si iv resonance lines in 
the more luminous stars, the Si iv result is based exclusively 
on the less luminous stars in our sample, and Vink et al. (2000) 
have already pointed out that there is a large discrepancy be- 
tween their predicted B star mass-loss rates and values based 
on Ha line profiles for stars with log(L/Lo)S 5.8. 

4. Implications for the nature of the mass-loss 

The major conclusions of the Letter are: 

1. The ionization structure of the winds of B supergiants is 
different from the O stars (Fig. 2). We caution, however, 
that most of these results are for low luminosity stars. 

2. There is no evidence of an abrupt change of the wind ion- 
ization in the temperature region of the bistability jump. 

3. The ionization fractions of non-CNO ions never exceed 
0.15 (for Al m) in the B supergiants - a result consistent 
with P v based results in the O stars (Massa et al. 2004). 

4. The ionization fractions of Si iv never exceeds 0.02. This 
result is based exclusively on lower luminosity B super- 
giants (where the Si iv lines are unsaturated) and holds 
when using either theoretical or empirical mass loss rates. 
It suggests that the theory is incomplete for lower luminos- 
ity stars and underscores the need for more empirical mass 
loss rates for these stars. 

For Galactic B supergiants, abundance anomalies cannot 
explain the very low {qi) for Al and Si, since these elements 
have well known abundances and are unaffected by stellar evo- 
lution. Revisions to the B star Teff scale are expected to have 
little affect on the predicted mass-loss rates at the temperatures 
where Si^"^, Si^"^ and Al^^ peak (see § 3). Note also that for this 
key temperature range, stellar radii based on the Humphreys 
& McElroy (1984) calibration differ from the Kurdtizki et al. 
(1999) radii by only ~ 5 - 15%. For the range of temperatures 
and ionizations considered in this Letter, it is highly unhkely 
that none of the ions are remotely close to dominant. As Massa 
et al. (2003) discuss concerning a similar result for P v in LMC 
O-type stars, we suspect that the small ion fractions indicate 



clumping in B supergiant winds. UV time-series spectroscopy 
has clearly demonstrated that the winds of early-B supergiants 
are highly variable and structured, with evidence for rotational 
modulation in many cases (e.g. Prinja et al. 2002), and opti- 
cal studies (see, Kaufer 1999) have reached the same conclu- 
sions for mid-to-late B supergiants. Furthermore, Blomme et al. 
(2002) find excess flux at millimetre wavelengths in e Ori (BO 
la) compared to homogeneous wind model predictions, indicat- 
ing considerable wind structure at several tens of stellar radii 
from the photosphere. Note, that if the winds are extremely 
clumped or structured, we must also re-evaluate the mass loss 
rates inferred from the radio and Ha observations. As a result, 
it is difficult at this time to infer how much of the low ioniza- 
tion fractions result from systematics in the measurement pro- 
cesses, and how much results from overestimates of the mass 
loss rates. 
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